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ABSTRACT 

We introduce the Making Galaxies in a Cosmological Context (MaGICC) program 
of smoothed particle hydrodynamics (SPH) simulations. We describe a parameter 
study of galaxy formation simulations of an L* galaxy that uses early stellar feedback 
combined with supernova feedback to match the stellar mass-halo mass relationship. 
While supernova feedback alone can reduce star formation enough to match the stellar 
mass-halo mass relationship, the galaxy forms too many stars before z = 2 to match 
the evolution seen using abundance matching. Our early stellar feedback is purely 
thermal and thus operates like a UV ionization source as well as providing some 
additional pressure from the radiation of massive, young stars. The early feedback 
heats gas to > 10 6 K before cooling to 10 4 K. The pressure from this hot gas creates a 
more extended disk and prevents more star formation prior to z = 1 than supernovae 
feedback alone. The resulting disk galaxy has a fiat rotation curve, an exponential 
surface brightness profile, and matches a wide range of disk scaling relationships. 
The disk forms from the inside-out with an increasing exponential scale length as the 
galaxy evolves. Overall, early stellar feedback helps to simulate galaxies that match 
observational results at low and high redshifts. 

Key words: galaxies: formation - galaxies: ISM - hydrodynamics - methods: N-body 
simulation 



1 INTRODUCTION 

As matter collapses in the early Universe, gas initially heats 
as kinetic energy is converted in to thermal dRees fc Ostrikerl 
1 19771 : iBirnboim fc Dekell l2003t) . In massive halos, the gas 
reaches temperatures where its cooling time becomes longer 
than the Hubble time (|Rees fc Ostrikei1ll977t ), which results 
in a theoretical maximum galaxy mass. However, in the cen- 
ters of galaxies, enough gas accumulates to increase the den- 
sity to the point that the gas efficiently radiates and cools. 
Such cooling is unstable since the radiation leaves through 
the optically thin surrounding hot gas. This process leads to 
overcooling since as the hot gas cools, it stops providing pres- 
sure support for the surrounding gas, more d ense gas is able 
to ac cumulate and cool even more efficiently (|White fc Reesl 
Il978ft . Th is problem is known as the "overcooling catastro- 
phe" (e.g. lBalogh et al.ll200ll ). 



Email: stinson 'at' mpia.de 



The "overcooling catastrophe" has long plagued 
simulations of disk galaxy formation. Fully cosmological 
numerical galaxy simulations consistently contain a mas- 
sive central concentration of stars (Navarro fc Bend 



199ll: iGovernato et al.l 
Scannapieco et al.l I2012F ) . 



20ld : 



12004 ; IStinson et aL 
This is evident in the cen- 
tral peak in rotation curves in simulated galaxies as well as 
in their surface brightness profiles. 

Further evidence for simulations forming too many stars 
comes from recent studies matching observations of galaxy 
stellar masses with their dark matter halo masses. One 
method is the abundance matc hing technique . To l ay the 
groundwork for this technique, IConrov et alj j2006) rank 
ordered halos by total halo mass from collisionless simu- 
lations and then did a rank ordering of galaxies based on 
their stellar mass from the Sloan Digital Sky Survey. The 
halos and gal a xies w ere then divided in several mass bins. 
IConrov et alj f2006) found that the correlation functions 
between halos in certain mass bins were well matched by 
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the correlation functions of galaxies in corresponding stel- 
lar mass bins. Several groups used this abundance match- 
ing technique to compare total halo masses with galaxy 



mg tecnmque to compare total naio masses witn galaxy 
stellar masses llConrov fe Wechslerll2009l ; lMoster et al'lboiol: 



IGuo et aill201oT Behroozi et alj2010h . Direct comparisons of 
stellar mass and total h alo mass have also been done based 
on s a tellite dynamics (iKlypin fe Pradal 120091; iMore et alj 
20091, 120111) and weak lensing |Mandelbaum et all 120091 ; 



Schulz et al]|2010! ). For a comparison o f all th ese techniques 
please see Figure 11 of lBehroozi et al" (l2010f). Each me thod 



shows a good level of correspondence. IGuo et alj (|2010h and 
ISawala et al.l (|201lT ) showed that nearly all the simulations of 
galaxy formation have formed many more stars than abun- 
dance matching predicts. 

Even without overcooling, Ivan den Bosch et all (|2002r ) 
showed that the amount of low angular momentum mate- 
rial in collapsed collisionless halos exceeds the amount of 
low angular momentum material observed in disk galaxies. 
Consequently, this low angular momentum material needs 
to be removed from the center of the system. 

Stellar feedback is the favor ed way of reducing star 
formation and laun c hing outflows (IScannapieco et aL 120081 : 
ISchave et all l2010h . iDekel fe Silkl ljl986h showed that su- 
pernova feedback can eject gas from galaxies with virial 
velocities up to 100 km s _1 . Semi-analytic models have 
found that significant amounts of stellar feedback are re- 
quired to match the low mass end of the luminosity 



quired to m atcn tne low mass end ot tne luminosity 
function (ISomerville fe Primackl Il999l; iBenson et al.1 120031: 



iBower et al.1 120061 . 120081 . l2012t ). ibutton fe van den Boschl 
l|2009h showed that stellar feedback can remove low angular 
momentum material. 

In hydrodynamical simulations, two methods are com- 
monly used to model stellar feedback. One is kinetic feed- 
back that adds velocity kicks to gas pa rticles to remove them 
from t he inner regions of galaxy disks dSpringel fe Hernquist 
2003; lOppenheimer fe Pavel 120061 : iDalla Vecchia fe Schavel 



20081 ). The other is thermal feedback in which stars sim- 



ply heat gas particles and allow the adia batic work of the 
parti c les to push other gas out of the way (iGerritsen fe 
19971: iThacker fe Couchmanl 120001 ; iKawata fe Gibson! 12003 
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2006). 



Since stars form in dense regions, the cooling times of 
the surrounding gas are short, and without help, th e gas 
will quickly radiate away all the supernova energy l|Kat j 
Il992l ). In real galaxies, the amount of gas necessary to exert 
a dynamical influence on the ISM small. In simulations, such 
small amounts of gas are difficult to model, so a common 
techniq ue has been to turn off cooling for a limited amount 
of time (IGerritsen fe IcQll997l;lThacker fe Couchmanll2000l : 



iBrook et al.ll2004 IStinson et alj|2006l ). 

It is not clear that kinetic and thermal feedback 
have significantly different effects than one another. 
iDurier fe Dalla Vecchial (|2012l ) showed that kinetic feedback 
has the same effect as thermal feedback when the hydrody- 
namics is left turned on. However, in many implementations 
of kinetic feedback, the hydrodynamic processes are disabled 
for a set period of time to maintain numerical converge nce 
|Springel fe Hernouisdl200"3l : lOppenheimer fe Davellioog) . 

It is also possible that thermal feedback does not need 
to rely on disabling cooling. IDalla Vecchia fe Schavel (|2012T ) 
recently showed that thermal feedback can have a significant 
dynamical effect if the total energy of all the supernovae ex- 



plosions is injected into one particle at one time. This large, 
one time energy deposition raises the particle's temperature 
high enough that its cooling time lengthens enough for the 
feedback to have a dynamical effect. 

For the first time in simulations. iGovernato et al.1 (|201dh 
showed that thermal stellar feedback can remove low angu- 
lar momentum material from dwarf gala xies, creating galax- 
ies with slowly rising rotation curves. IBrook et all (|201lD 
showed that the ejected low angular momentum material 
becomes part of a galactic fountain and can be reaccreted 
onto the disk with higher angular momentum than it left. 
ISawala et ail i|201ll ) showed that even using a high density 
threshold for star formation that made the stellar feedback 
more efficient, these simulations produced more stars than 
is predicted by the stellar mass-halo mass relationship. 

Three other recent simulatio ns have shown succ ess at 
forming realistic L* disk galaxies. iGuedes et al.l l|201ll ) used 
high resolution and a high amou nt of supernova feed back 
to produce a realistic disk galaxy. Ag ertz et al.l l|201ll ) sug- 
gested that using a low star forming efficiency is helpful for 
creating extended disks. They also use models that increase 
the amount of supernova energy deposited above the canon- 
ical 10 51 erg per supernova. These simulations do the best 
job flattening out the central region of their rotation curve. 
iMcCarthv et al.1 (|2012l ) also form galaxies with flat rotation 
curves in a large, low resolution cosmological volume. They 
employ a kinetic feedback scheme that uses les s than the 
canon ical 10 51 erg supernova energy, though see iKav et ail 
(2002) for a demonstration of how kinetic feedback provides 
a larger effect than thermal feedback. The galaxies each con- 
tain a few thousand particles and employ a gravitational 
softening of ~ 1 kpc, which makes an examination of the 
detailed structure of the galaxies challenging. 

These simulations show that it is possible that more 
feedback than the canonical 10 51 erg is required. Instead of 
simply increasing the amount of energy released b y super- 
novae feedback, we note that iMurrav et al.l l|2010l ) showed 
that there can be significant feed back effects from stars be- 
fore they explode as supernovae. iHopkins et al.l l |201ll ) in- 
corporated a kinetic radiation pressure feedback into sim- 
ulations of isolated disk galaxies and found that the feed- 
back could strongly regulate star formation. While we do 
not have the resolution in cosmological simulations to im- 
plement a similar kinetic scheme, we implement a scheme 
based on thermal pressure to provide feedback during the 
time between when stars are formed and the first SN star 
exploding. 

Using this feedback prescription, iMaccio et all (|2012t) 
showed that the feedback removes low angular momentum 
dark matter in galaxies up t o nearly L*, p r oduci ng cored 
dark matter density profiles. IStinson et al.l (|2012f l showed 
that the metal rich outflows created by this feedback match 
observations of O VI in the cir cum-galactic medium of star 
forming galaxies. IBrook et al.l (2012b) also showed that a 
sample of lower mass galaxies form disks that follow a wide 
range of disk scaling relationships. 

Here, we present a detailed study of how varying the 
key free parameters in our simulations affect the morphology 
and evolution of galaxies. In Sj2] we describe the galaxy we 
model and the physics used in the simulation. In JJ] we show 
how varying parameters affects the mass of stars formed, the 
morphology of the galaxy, and how the galaxy evolved. 
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2 SIMULATIONS 

We use gl536, a cosmological zoom simulations drawn from 
the McMaster Unbia sed Galaxy Simulatio ns (MUGS), in our 
parameter study. See lStinson et all (|2010t ) for a complete de- 
scription of the creation of the initial conditions. gl536 has 
a total mass of 7 x 10 11 Mq, a spin parameter of 0.017, 
and a last major merger at z — 2.9. Using the physics em- 
ployed in the original MUGS simulations, gl536 also formed 
an exponential disk with an classical bulge with a central 
surface brightness of \ii — 14, a total face-on magnitude of 
M r =-21.7 and g - r colour of 0.62. 

Using the smoothed part i cle h ydrodynamics (SPH) 
code GASOLINE l|Wadslev et al.l [20o3 ). we varied the stel- 
lar feedba ck to reduce s tar fo rmation to the amount pre- 
scribed in iMoster et all |2012T ). We first increased the en- 
ergy input from supernovae to 10 51 erg and changed the 
in itial stellar ma s s func tion (IMF) from what was described 
in lKroupa et alj (|l993l ) t o an IMF t hat co ntains more mas- 
sive stars as described in IChabrierj (|2003h . Neither of these 
changes sufficiently limited star formation, so we turned to a 
source of energy that is often ignored in cosmological simu- 
lations, the luminosity emitted by massive s tars before they 
explode as supernovae. iHopkins etafl (|201lh has shown that 
the inclusion of radiation pressure in high resolution simula- 
tions can reduce star formation. The supernovae explosions 
in our previous simulations did not happen until 4 Myr af- 
ter the formation of a stellar population. However, molecular 
clouds are disrupted much sooner than this, so the inclusion 
of this early stellar feedback is necessary to reduce star for- 
mation before supernovae start exploding. 

2.1 Gas Cooling 

In addition to energy input, the gas cooling has a large effect 
on how many stars fo rm. The cooling u sed in this paper is 
described in detail in IShen et alj J2010 | ). It was c alculated 
using CLOUDY (version 07.02; iFerland et al.1 Jl998l)') includ- 
ing ph otoionization and heating from the lHaardt fc Madaul 
|2005T ) UV background, Compton cooling, and hydrogen, he - 
lium and metal cooling from 10 to 10 9 K. IShen et all ifcOlOh 
showed that T > 10 4 K metal cooling including the pho- 
toionization effects of the UV background reduces cooling 
by a n order of magnitude from purely collisional cooling 
rates ISutherland fe Dopital (|l993l l. 

In the dense, interstellar medium gas, we do not impose 
any shielding from the extragalactic UV field. While the 
cosmic UV background used is lower than the typical UV 
field within galaxies, it does provide heating similar to that 
expected in the ISM. 

When metal cooling is included below 10 4 K, it leads to 
the well-known tw o-phase instability of the ISM (Field 1965; 
iField et al.1 1 19691} . which produces densities over 100 cm 
and temperatures of 100 K in these simulations. These values 
correspond to a Jeans mass of <6000 Mq and a Jeans length 
of <8 pc. Each of these values are far below both the gas 
particle mass, 2x 10 5 Mq, or 310 pc softening length. In real- 
ity, this phase would become even denser, self-shielded, and 
ultimately molecular, but we cannot resolve this process and 
must therefore estimate the star formation rate from the re- 
solved dense gas. To avoid issues with Jeans fragmentation, 
a minimum pressure is established in the gas as described in 



iRobertson fc Kravtsovl ((2008) . A maximum density is also 
set to 400 cm~ 3 using a minimum smoothing length ~ of the 
310 pc gravitational softening. 

2.2 Star Formation and Feedback 

The simulat ions use a common recipe for star formation 
described in IStinson et al.l (|2006r ) that we summarize here. 
Stars form from cool (T < 15, 000 K), dense gas. The metal 
cooling readily produces dense gas, so the star formation 
density threshold is set to the maximum density at which 
gravitational instabilities can be resolved, 32 j""* (nth > 9-3 
cm -3 ), where M gas = 2.2 x 10 5 Mq and e is the gravitational 
softening (310 pc). Such gas is converted to stars according 
to the equation 

AM* ={ _ M gas 

At tdyn 

Here, AM* is the mass of the star particle formed, At is 
the timestep between star formation events, 8 x 10 5 yr in 
these simulations, and tdyn is the gas particle's dynamical 
time, c* is the efficiency of star formation, in other words, 
the fraction of gas that will be converted into stars during 

tdyn - 

Stars feed both energy and metals back into the in- 
terstellar medium gas surrounding the region where they 
formed. Supernova feedback i s implemented using the blast- 
wave formalism described in IStinson et al.1 (|2006T l. In this, 
Type II supernovae are assumed to explode due to the 
core collapse at the end of lifetime of stars greater than 
8 Mq. Stellar lifeti mes are based on the Padua lifetimes 
|Raiteri et al.| [l996). Since the gas receiving the energy is 
dense, it would quickly be radiated away due to its efficient 
cooling. For this reason, cooling is disabled for particles in- 
side the blast region 

in 1 - 74 rr0.32 -0.16 5-0.20 / n \ 

rcso = 10 E 51 n P 04 pc (2) 
and for the length of time 

. i n 6.85 rr.0.32 0.34 5-0.70 / Q \ 

tcso = 10 £ 5 i n P 04 yr (3) 
given in iMcKee fc Ostrikeij (|l977r i. 

2.3 Early Stellar Feedback 

In an effort to provide a more physical model for stellar feed- 
back, radiation energy from massive stars is now considered. 
Heating is introduced immediately after stars form based 
on how much star light is radiated. The early energy input 
from high mass stars has been entirely ignored in most pre- 
vious cosmological simulations. In such cases, 4 Myr pass 
after stars form before the first supernov a explodes (some 
simulations use even longer delays, see ISlvz et all 120051 : 
iDalla Vecchia fc Sch avc 2012), during which time stars can 
con tinue forming w i thout any affect of stellar feedback. 

IHopkins et al.1 |201ll ') describes how such energy could 
be considered using a kinetic scheme where the radiation 
pressure drives winds out of massive star clusters. Our sim- 
ulations do not have the resolution to use kinetic feedback 
because it requires the choice of a direction. A kinetic im- 
plementation of feedback requires many particles to repre- 
sent a molecular cloud so that it can be torn apart. In our 
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Table X. Simulation data 



Name 


color 


c* a 




Ctd C 


IMF d 


M* [M ] e 


M v f 


Fiducial 


red 


0.1 


0.1 


0.05 


C 


2.3 x 10 10 


-20.6 


MUGS 


black 


0.05 





0.05 9 


K 


8.3 x 10 10 


-21.4 


Low Diffusion 


green 


0.1 


0.1 


0.019 


C 


3.0 x 10 10 


-21.1 


High Diffusion 


yellow 


0.1 


0.175 


0.05 s 


c 


2.5 x 10 10 


-21.4 


No ESF 


magenta 


0.1 





0.05 


c 


3.4 x 10 10 


-19.5 


120% SN energy 


cyan 


0.1 





0.05 


c 


1.8 x 10 10 


-20.0 


High ESF 


blue 


0.1 


0.125 


0.05 


c 


1.1 x 10 10 


-19.9 



a Star forming efficiency 

k Early Stellar Feedback Efficiency 

c Thermal Diffusion coefficient 

d Initial Mass Function: C= IChabrierl l2003f) ; K= lKroupa et al.l jl993h 

e M, is the total stellar mass 
f V-band magnitude 

9 All particles diffused thermal energy including those with cooling shut off 



simulations, only one or a few particles represent a molecu- 
lar cloud, so we utilize thermal pressure, which is isotropic. 
Using thermal feedback provides pressure support and in- 
creases gas temperatures above the star formation threshold 
to decrease star formation. We do not disable cooling during 
these early times since massive young stars radiate a large 
amount of energy. Thus, after heating the gas to T > 10 6 K, 
the gas rapidly cools to 10 4 K, which creates a lower den- 
sity medium than if the gas were allowed to continue cooling 
until supernovae exploded. 

To model the luminosity of stars, a simple fit of the 
mass-lu minosity relationshi p observed in binary star sys- 
tems bv lTorres et"all i|2010l ) is used: 

_L J to)' 2 > M < WM ° (4) 

Lq I 100 (l^) 2 ,M>WM e 

Typically, this relationship leads to 2 x 10 50 erg of energy 
being released from the high mass stars per Mq of the entire 
stellar population over the 4 Myr between the star's forma- 
tion and the commencing of SNII (SNII inject ~ 10 49 erg 

per M a ). 

iLeitherer etafl (l999) found that approximately 10% of 
the total stellar flux is emitted in the UV, so we use with 
an ea rly stellar feedback efficiency, e es /=10%. l|Frever et al.l 
2006) found that stellar photons do not couple efficiently 
with the surrounding ISM. Correspondingly, radiative cool- 
ing is not turned off for this form of energy input. The dense 
gas thus cools rapidly to 10 4 K in the ISM, which mimics 
the formation an HII region. Though the dynamical effect 
is minimal, early stellar feedback effectively halts star for- 
mation in the region immediately surrounding a recently 
formed star. Using e es f = 10% limits star formation to the 
amount prescribed by the stellar mass-halo mass relation- 
ship. 

2.4 Diffusion 

lAgertz et al.l ([2007) highlighted how long cold gas clumps 
can survive in a surrounding bath of hot gas in SPH rela- 
tive to Eulerian grid simulations. In an effort to decrease 



this mixing time, IWadslev et al.l l|200cf ) implemented dif- 
fusion into SPH based on mixing that should result in 
shearing flows. The amount of mixing depends on the 
magnitude of the local velocity shear field multiplied by 
the square of a length compara ble to the grid or resolu- 
tion scale. IWadslev et al.l (|2008l ) showed that diffusion of 
ther mal energy e f fectiv ely limits the lifetime of the blobs 
thatlAeertz et alJ <|2007t > described. In MUGS ijStinson et all 
120101 ). both metals and heat diffused using this scheme with 
a diffusion coefficient of 0.05. The most striking effect of 
the thermal diffusion was that cold accretion filaments were 
heated and expanded. Cold blobs persisted in these simula- 
tions. Results of varying the thermal diffusion coefficient are 
presented in §3] 

During this study, it was determined that the heat diffu- 
sion severely reduced the efficiency of the adiabatic feedback 
scheme to drive outflows. Hot particles with their cooling 
shut off frequently passed cooler particles with their cooling 
turned on. Heat would diffuse from the hot to cold parti- 
cles, reducing the temperature of the particles and thus the 
efficiency of the feedback. This would happen mostly as hot 
particles orbited around the disk and occasionally after par- 
ticles had been ejected from the disk. Since the feedback 
is dependent on particles not cooling, we disabled thermal 
diffusion in interactions between particles in which either 
particle had its cooling shut off. We note that disabled cool- 
ing has the biggest effect in the dense disk. We find less than 
1% of hot (T> 10 4 K) particles with their cooling disabled 
are outside the disk at z = and such particles are never 
found beyond 25 kpc at any time. 

Figure [1] shows how the star formation histories of the 
different implementations of thermal diffusion vary. In the 
fiducial model, the star formation declines until the end of 
the simulation. The scheme where all gas particles always 
diffused required a 75% increase in cesf to limit the star 
formation to the stellar mass-halo mass relationship. The 
shape of the star formation history for the high diffusion 
model is different. It increases throughout the history of the 
galaxy and begins to start a starburst over the last 1 Gyr. 
This shows how with too much thermal diffusion star for- 
mation acts as a positive feedback: the more stars form, the 
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Fiducial 

Low Diffusion 

High Diffusion 




11 12 

Time [Gyr] 

Figure 1. The star formation history of galaxies simulated with 
varying amounts of thermal diffusion. Here, the simulation evolves 
from left to right. This plot focuses on the last 5 Gyr of the simu- 
lations since that is when the biggest difference in star formation 
rate is seen. 



more gas cools enabling more star formation. The low ther- 
mal diffusion simulation (green) demonstrates a milder in- 
crease in star formation rate at the end of the simulation, but 
it remains above the fiducial value. An examination of the 
simulation reveals that this is due to cold gas clouds c ooling 
out of the hot halo as descr ibed in lMaller fc Bullock! (|2004l l 
and iKaufmann et alj |2009h . In the revised treatment, the 
limited diffusion helps to limit the formation of cold blobs 
and reduces the star formation rate at late times in the sim- 
ulation. 



3 RESULTS 

To find an optimal star formation and feedback recipe, we 
attempt to match the stellar mass-halo mass relationship. 
The parameters for all of our simulations are presented in 
Table 1. Figure [2] shows a mock observation of the fiducial 
galaxy. The image is 50 kpc on a side and was creat ed using 
the M onte Carlo radiative transfer code sunrise (|jonssonl 
2006). The image bri ghtness and contrast are scaled using 
asinh as described in iLupton et al. (2004) since disks have 
an exponential surface brightness profiles meaning the im- 
ages span a wide range in surface brightness. The image 
shows a thin, young stellar disk with clumpy and extended 
star formation, a dust lane, and a red bulge component 



3.1 Total Mass 

Figure shows the mass of the stellar component at z — 
as a function of halo mass. For the sake of clarity we 
have limited the number of simulations displayed to those 
which demonstrate our walk through parameter space. The 
coloured stars show the z — stellar and halo masses of 



a Movies of the simulations can be found 
http://www.mpia.de/~stinson/magicc 



at 



Figure 2. Face-on and edge-on images of the fiducial galaxy 
at z = 0. The images are 50 kpc on a side and were created 
using the Monte Carlo radiative transfer code SUNRISE. The image 
brightness and cont rast are scaled using asinh as described in 
ILupton et all l|2004h . 



the simulations of gT536 used in this parameter study. Each 
2 = point is accompanied by a line that shows the path 
along which it evolved through this diagram. A couple of 
lines are provided for reference. The dotted line is the cos- 
mic baryon fraction of the total halo mass. The dashed line is 
the ste llar mass halo mass relationship found in lMoster et al.l 
|2012[ ). The black plus signs are the old results from MUGS. 
These lower feedback runs are similar to many other low 
feedback simulations where galaxies turn about half of a 
galaxy's baryons, as predicted by the cosmic baryon frac- 
tion, into stars. 

The red star represents the best fit to the stellar mass- 
halo mass relationship. It includes early stellar feedback with 
thermal diffusion disabled for particles with their cooling 
disabled. It uses a star formation efficiency, c* =0.1, an early 
stellar feedback efficiency, e es / =0.1, and a thermal diffusion 
coefficient, C t d =0.05. The fiducial simulation ev o lves a long 
a path generally below the z = iMoster et alj |2012T ) re- 
lationship, which is shown in i|3.3l to be similar to how the 
observed stellar to halo mass ratio evolves. A small change 
in the amount of feedback causes a significant change in the 
number of stars that form. When e es / increases from 0.1 to 
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Figure 3. Stellar mass plott ed as a function of halo mass. The 
abundance matching fit from Most er et al. I <2012h and the cos- 
mic baryon fraction lines are shown as reference. All the MUGS 
simulations formed too many stars along a line where half the 
baryons have turned into stars and half are in the hot gas halo. 
The simulations of gl536 are shown as coloured stars and lines. 
The legend provides a brief description of which simulation each 
of these colours represents, while a more extensive description can 
be found in Table 1. The fiducia l simu lation is coloured red and 
lies closest to the iMoster et al. relationship. The yellow 

star representing the simulation where all gas particles diffuse 
thermal energy also lies close to the relationship line, but it is 
rapidly forming stars and will quickly overshoot the relationship. 
The blue point uses the same physics as the fiducial run, but in- 
creases the e e3 f by 25% and forms far fewer stars. The magenta 
line shows a simulation run without early stellar feedback. While 
it ends up close to the fiducial run at z = 0, it follows a dramati- 
cally different evolutionary path. 

0.125 (represented in blue), the stellar mass decreases by a 
factor of 2. 

The magenta star and line represent the simulation in 
which only the supernova feedback was increased from the 
MUGS value of 4 x 10 50 erg per supernova to 10 51 erg. In 
this simulation, the thermal diffusion is turned off for par- 
ticles with their cooling disabled. This change significantly 
decreases the mass of stars forme d , thou gh the stellar mass 
still lies above the IMoster et aD <|2012f > relationship. This 
motivated us to add early stellar feedback. 

As a check to see whether the addition of more energy 
or the timing of that energy addition had a stronger effect, 
we ran a series of simulations without the 10 51 erg per su- 
pernova energy constraint. The best result from that series 
is represented by the cyan line and star. It shows that a 
20% increas e of the supernova energy lowered the stellar 
mass to the IMoster et al.l (|2012 | ) relationsh i p. Its evolution 
more closely follows the z = IMoster et all (|2012l ) relation- 
ship than the simulation using 10 s1 erg. A gain, we refer the 
reader to ^3.3lto see that the evolu tion does not follow what 
is described in IMoster et al] (|2012h . 

The yellow star and line show the reduced effect of stel- 
lar feedback when thermal diffusion is active for particles 
with their cooling shutoff. This requires the early stellar 
feedback efficiency, e es /, to be increased to 0.175 to produce 
a galaxy that matches the stellar mass-halo mass relation- 



Figure 4. B-band surface brightness profile for gl536 simulated 
at a variety of feedback strengths. Each surface brightness profile 
has an exponential component with a central surface brightness 
around 21 mag arcsec -2 . All the simulations that contain too 
many stars in Figure \3\ have an exceptionally steep central pro- 
file that ends at a central surface brightness brighter than most 
observed galaxies. The galaxy that uses the MaGICC feedback 
deviates slightly from the pure exponential line, but does not in- 
crease nearly as quickly as the simulations with lower feedback. 

ship at z = 0. However, its evolution is much different from 
the other galaxies. Its stellar mass is increasing sharply at 
2 = 0. Figure [1] shows that its star formation rate is also 
dramatically increasing, so it will quickly overshoot the re- 
lationship. The drastic increase in star formation rate is the 
result of thermal feedback diffusing away instead of con- 
tributing to launching winds. 

Simply reducing the thermal diffusion for all gas parti- 
cles from 0.05 to 0.01, as shown by the green star, required 
a corresponding reduction to e es f to form enough stars to 
match the stellar mass-halo mass relationship. However, the 
reduced thermal diffusion allowed cold clouds to condense 
out of the hot halo at the end of the simulation, which again 
led to a increasing star formation rate, though not as dra- 
matic as in the higher thermal diffusion case. 

As we made our parameter study, we found a degener- 
acy between c* and e es f, so that a simulation where c*=0.05 
and e es /=0.15 produces the same results as a simulation 
where c*=0.1 and e es /=0.1. Such a degeneracy is not sur- 
prising as the early feedback is immediate and limits star 
formation efficiency. 

3.2 Mass Distribution 

Having used the stellar mass-halo mass relationship as 
our constraint, we now examine other properties of the 
galaxies. Figure [4] shows surface brightness profiles pro- 
duced using different amounts of stellar feedback in gl536. 
The one dimensional, face-on surface brightness profiles 
were constructed us ing the stellar population models from 
iGirardi et al.l |2010). Only stars in a disk with a radius of 
15 kpc and height 2 kpc above and below the disk midplane 
are included. The profiles represent the mean values of 100 
linearly spaced azimuthal bins. An exponential component 
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Figure 5. The evolution of the B-band surface brightness profile 
in the fiducial simulation. 




r[kpc] 



Figure 6. Plot of v c = \J GM(r)/r as a function of r for the all 
the models of gl536 and g5664 at z = modeled with varying 
amounts of stellar feedback. 



of the disk was fit using stars with radii greater than 5 kpc. 
These fits are represented by the dashed lines the same color 
as the surface brightness profiles. Each exponential fit has a 
scale length around 4 kpc and central surface bri ghtnesses 
aroun d 21 in the B-band, slightly brighter than the lFreemanl 
i|l970l ) Law. 

The fiducial simulation forms the right amount of stars 
and has an exponential surface brightness profile that steep- 
ens at 5 kpc to reach a central B surface brightness of 
17 mag arcsec -2 . By contrast, the simulations of gl536 
without early stellar feedback each have exponential sur- 
face brightness profiles with slightly longer scale lengths 
but more extensive bulges. These bulges are characterized 
by higher Sersic index components in their inner regions 
that reach brighter central surface brightnesses than the 
simulations with early stellar feedback. This high central 
surface brightness is a common symptom of excess stellar 
mass in the central galactic region found in many galaxy 
simulations (S carmapie co et alj l201fJ : IStinson et al.1 bold ; 
IScannapieco et al]|2012ft . 

Figure[S]shows the evolution of the B surface brightness 
profile in our fiducial model. It shows the disk growing from 
the inside out. For the bulk of its evolution, the surface 
brightness profile remains exponential. In the final two time 
steps shown, an excess develops in the center. The scale 
length of the disk undergoes an evolution from 1 kpc at 
z — 2 to 4 kpc at z = 0. 

To examine the underlying mass distribution of the 
galaxy, Figure [6] shows the rotation curves, the circular ve- 
locity, v c — \J GM(r)/r, as a function of radius, for the same 
galaxy models. The fiducial model has a nearly flat rotation 
curve. The rotation curve for the MUGS gl536 has a large 
central peak. The simulation without early stellar feedback 
also exhibits a high central peak in its rotation curve. How- 
ever, the simulation with 20% more supernova energy has a 
much flatter rotation curve, though not as flat as the fiducial 
model. 

The simulations with too much feedback have slowly 
rising rotation curves that are in conflict with the nearly 
flat rotation curve observed for the Milky Way and other 




star 

°0 2 4 6 8 10 12 14 16 

r / kpc 

Figure 7. Plot of v c = ^jGM(r)/r as a function of r for the 
simulated galaxy at z = 0. Each component, dark matter (solid 
red), gas (green, dotted), and stars (light blue, dashed), is plotted 
separately to show the matter distributions of each component. 

L* galaxies (jCourtead 119971 ; Ide Blok et al.ll2008l ). We have 
drawn sample model fits for the fi ducial and high ear ly stel- 
lar feedback simulations following iReves et all (|201ll ). They 
use Rto values, radius at which the profile turns over, of 0.25 
and 1 kpc, respectively. The slowly rising rotation curves 
sh ow the best agreement with the arctangent models used 
in lReves et al.1 (|201ll ). where typical Rto values are approx- 
imately the same as the disk scale length for galaxies with 
stellar masses similar to gl536's 2 x 10 10 M , which means 
that the rotation curves do not reach their peak until well 
beyond the disk scale length. 

Even though the rotation curves slowly rise in the plot- 
ted region, the peak of the rotation curve, 180 km s _1 is still 
higher than the value of v V i r , 110 km s _1 . In other words, 
the rotation curves turn over and decline outside the plotted 
region. 
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Figure 8. Peak circ ular velocity as a fu nction of r-band magni- 
tude compared to the lReves et al,l l l201lf ) observations. The simu- 
lations are plotted with the standard colour scheme used through- 
out this paper. The fiducial simulation (red) lies directly in the 
center of the observations. The simulations without early stellar 
feedback lie significantly above the observations. 
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Figure 9. The star formation history of galaxies simulated with 
varying amounts of stellar feedback. Here, the simulation starts at 
Time = Gyr, on the left of this plot. The star formation rate is 
drastically reduced with the increased stellar feedback. The shape 
also changes from a quick rise and exponential falloff to a flatter 
evolution. 



The rotation curve for the fiducial galaxy is separated 
into the contributions from its constituent particle types in 
Figure [7] Stars dominate the potential in the central regions 
of the galaxy. The curves represent the mass included from 
just that component, so it does not represent the speed at 
which those particles are orbiting. Each component will orbit 
close to the total velocity. 

The rotation velocity information enables another com- 
parison between the amount of stars that have formed 
and the mass of the galaxy, the Tully-Fisher relationship 
jTullv fc Fished 1 19771 '). The Tully-Fisher relationship com- 
pares the luminosity of a galaxy with its circular velocity. 
It is possible to argue for the choice of many different radii 
at wh ich to measure the circular velocity. In IStinson et al.l 
(2010), the galaxies were reported as fitting the Tully-Fisher 
relationship because the circular velocities were taken at 2.2 
Rd, where Rd is the scale radius, which was outside the cen- 
tral velocity pea k. Here, we instead choose to compare our 
simulations with lReves et al] (|201ll ). who use an arctangent 
model and report the peak velocities found using that fit. We 
compare these with the peak velocity from our simulations. 

The fiducial simulation along with the other simula- 
tions that use early stellar feedback lie in the middle of the 
observed relationship. One exception is the high diffusion 
simulation, which has a bright r-band magnitude due to 
its significant late star formation but still has a low peak 
velocity. The simulations without early stellar feedback lie 
above most of the galaxies in the observed relationship due 
to their high central velocity peaks. No ESF and MUGS lie 
well above the observations, but 120% SN energy lies only 
slightly above the observations. 



3.3 Star Formation History 

The star formation history in Figure [9] shows that vary- 
ing feedback has a significant impact on when stars form. 



The MUGS case shows a fairly standard star forma- 
tion history for most simula tions reported in the litera- 
ture l|Scannapieco et aU l2012h . Star formation follows the 
buildup of total halo mass as shown in Figure 1111 All of 
the simulations without early stellar feedback experience a 
fairly broad peak of star formation just after 2 Gyr. This 
is the time period when many nearly equal mass mergers 
are rapidly building the halo mass. When the supernova en- 
ergy is increased to 120%, star formation after this peak is 
kept below 3 Mq yr" 1 . Introducing the early stellar feed- 
back reduces or eliminates this burst from the simulations 
shown here because the early stellar feedback prevents gas 
from collapsing into the central 2 kpc as shown in S|4] Fol- 
lowing the early peaks of star formation the star formation 
rate exponentially declines with intermittent bursts during 
mergers. 

The star formation rate has been calculated using 100 
bins each 137 Myr wide in Figure [9] Figure [T] shows the star 
formation rate calculated using 25 Myr bins, which shows 
the bursty nature of star formation that is characteristic 
of all of the galaxies outside of the MUGS. There are mi- 
nor episodic bursts every 200-300 Myr overlaid onto the un- 
derlying shape of the star formation histories. These minor 
bursts correspond to the dynamical time of the galaxy, ap- 
proximately the time it takes for gas to collapse after it is 
heated by stellar feedback (this phenomenon is descr ibed in 
more detail for dwarf galaxies in IStinson et afll2007h . 

While the ESF galaxy does not show an increase in star 
formation at z — 4 (2 Gyr), all the galaxies show a peak of 
star formation at z — 2 (3.75 Gyr). This corresponds to the 
time of the first merger of halos with masses > 10 11 Mq. 
The "no ESF" simulation (magenta line) with Esn = 10 51 
erg but without early stellar feedback undergoes a moder- 
ately long, 2 Gyr starburst due to this merger. The fiducial 
simulation with early stellar feedback (red) has a shorter 
starburst (~ 100 Myr) and smaller amplitude (8 Mq yr -1 ) 
at this time. We note that a large population of small emis- 



© 0000 RAS, MNRAS 000, 000-000 



MaGICC Parameter Study 9 



14 3 



0.04 



0.03 



0.02 



0.01 



Moster+ (2012 
Fiducial 

120% SN energy 
Low Diffusion 
High Diffusion 
No ESF 
High ESF 




time [Gyr] 

Figure 10. The evolution of the ratio of stellar mass to total 
halo mass for simulations using varying amounts of stellar feed- 
back. Here, the simulation star ts at Time = Gyr , on the left of 
this plot. The black line is from lMoster et al.l j2012ll that includes 
a detailed treatment of substructure to make the relationship be- 
tween halo and stellar mass. The grey region represents the lcr 
variation in the model. The simulations evolve mostly within the 
lcr variation. 



sion line galaxies has been detected at z = 2 that have a 
major starbu rst as their first sign of significant star forma- 
tion activity (|van der Wei et all 20 111 ). 

Shortly after the major merger, the fiducial simulation 
has another peak at 4 Gyr that corresponds to the accretion 
of gas remaining from the z — 2 merger. Its peak is lower 
and it lasts for a longer time. There are a couple of small, 
gas rich mergers that create a second star formation rate 
peak at 7 Gyr. This ends the merging history for gl536 and 
for the most part, the star formation rate slowly declines as 
star formation uses up the gas in the disk. The simulations 
without early stellar feedback decline more rapidly than the 
fiducial simulation and have a lower star formation rate for 
the last 7 Gyr of the simulation. This is due to those galax- 
ies either using up their gas reservoirs faster (No ESF) or 
ejecting the gas further so that it is no longer available for 
star formation (120% SN energy). 

The star formation history of the fiducial simula- 
tion compares well to the star formation history of the 
Milky Way, though such studies have necessarily fo- 
cused on the solar neighborhoo d llRocha-Pinto fc Maciell 



19971; iRocha-Pinto et all |2000| ; iHernandez et all I2006T : 



Cignoni et al. 20061 ). These studies find that the solar neigh- 



borhood has a mostly flat star formation rate with factors of 
a few var iation and a peak a t 2 ~ 1. Looking at nearby dwarf 
galaxies, IWeisz et all l|201ll ) finds that star formation histo- 
ries are also mostly flat with variations limited to factors of 
a few, which again compare s well with the sim ulation. We 
note that the galaxy in the Wei sz et al.l l|201ll ) sample are 
all at least 3 magnitudes fainter than the simulated galax- 
ies and have star formation rates two orders of magnitude 
lower. So, in a general sense, the star formation history of 
the fiducial simulation agrees well with the observations. 

Figure [10] presents another view of the evolution of stel- 
lar mass. It shows the evolution of the ratio of stellar mass 



to total halo mass for simulations using varying amounts of 
stellar feedback. The black line shows the evolution of the 
ratio based on the evolution of the total halo mass inside r2oo 
according t he recently publishe d stellar mass-halo m ass rela- 
tionship in lMoster et all (|2012l ). lMoster et all (|2012l ) adds a 
detailed treatment of substructure to match the abundances 
between dark matter halo and stellar mass as well as using 
updated high redshift luminosity functions. The black line is 
calculated using the mass evolution of the gl536 halo includ- 
ing baryons. The grey region represents the lcr variation for 
the model. The simulations that use early stellar feedback 
evolve mostly within the lcr variation. The fiducial (red) 
simulation follows the black line most closely throughout 
the evolution of the galaxy. The simulations without early 
stellar feedback evolve outside the lcr variation early in their 
evolution because of the large amount of early star forma- 
tion. A quick glance at Figure [9] indicates that the MUGS 
simulation that follows a star formation history typical of 
many published simulations would range even further out- 
side the acceptable range of stellar masses at early times. 
The 120% SN energy simulation does not evolve very far 
outside the variance region, but it also does not follow the 
evolution of the mean line closely at all. It is well above the 
mean for the first 7 Gyr and well below the mean for the 
last 5 Gyr. It is possible that this unique evolution has to do 
with the uniquely quiet merger history of gl536, but it could 
also indicate a mismatch. The high ESF model also evolves 
outside of the observed range of galaxy stellar m asses. It 
forms fewer stars at z = than lMoster et ail (|2012l ) predict. 

3.4 Mass Evolution 

It is also possible to track the mass history for the other 
components of the galaxy. Figure QT] shows the mass history 
for four of the simulations. For each, the blue line represents 
the total galaxy mass inside r V i r , where r„; r is the radius at 
which ~fh^a = 390pback g round- The black line is at one-sixth of 
the total mass, representing the cosmic fraction of baryons. 
The magenta line shows the actual mass of baryons inside 
the galaxy. The dashed cyan line represents the mass of gas 
with T< 10 5 K while the dash-dotted red line represents 
the mass for T> 10 5 K, and the dotted green line shows the 
evolution of stellar mass. The total mass line clearly shows 
the two significant mergers at z = 2 and 2 = 1. 

The origin of hot gas mass is coincident with the on- 
set of star formation, indicating that stellar feedback could 
play a role in creating the hot halo; it may n ot o nly be accre- 
tion shocks as described in iToft et all l|2002i ) and[c rain et al.l 
(2010). Proving this claim is beyond the scope of this paper, 
so we leave it for future work. Both stellar mass and hot gas 
mass increase significantly during the z — 2 merger, though 
the star formation increase is greatest in the case without 
early stellar feedback. The no ESF simulation also evolves 
to an equipartition between the three baryonic phases, stars, 
hot gas, and cold gas. It appears that the other simulations 
are converging to this equipartition, but it is hard to predict 
whether the stars will reach the same value as hot and cold 
gas. 

In the two upper panels of simulations without ther- 
mal diffusion in particles with their cooling shut off, 25% 
of the baryons escape the potential well of the halo during 
the 2 = 2 merger. Some baryons also initially escape in the 



© 0000 RAS, MNRAS 000, 000-000 



10 Stinson et al. 



513 2 1 



513 2 1 



a 10 1 



10 1 



a 10 1 



10 1 









/" Fiducial 




* High ESF 






" - - --,-.■=.-»—--- — - 




/ 

/ r ' 




Hi / 






✓ High Diffusion 




✓ No ESF 












Total Mass - ■ - 






/ 6 =0.16 _ _ M^^a 

Stellar Mass ^ 



2 4 6 8 10 12 
time [Gyr] 



2 4 6 8 10 12 
time [Gyr] 



Figure 11. The mass evolution of the total halo mass (blue), 
stellar mass (dotted- green), cold gas (T< 10 s K, dashed cyan), 
and hot (T> 10 5 K, dash-dotted red) gas within r v i r for the 4 
simulations. The thin black line shows the cosmic baryon fraction 
of the total halo mass, while the thin magenta line shows the 
actual baryonic mass content of the halo. 



simulations with thermal diffusion between all particles, but 
then baryons reaccrete. 



4 DISCUSSION 

The question remains about how the early stellar feed- 
back more successfully reproduces observed galaxy proper- 
ties than models without early stellar feedback. The greatest 
difference seems to be in the galaxy centers; the simulations 
with early stellar feedback have flat rotation curves while 
those without have central peaks. To examine this, Figure 
1121 compares the evolution of the inner 2 kpc of the fiducial 
simulation with the 120% SN energy simulation. In Figure 
it was shown that the 120% SN energy simulation forms 
more stars early throughout the galaxy than the fiducial 
simulation. Figure [12] shows that this pattern is replicated 
in the central 2 kpc, where most of the star formation is oc- 
curring. Figure [12] also shows that the higher star formation 
follows the higher mass of cold gas in the central region of 
the 120% SN energy simulation. So, there is a clear differ- 
ence in the cold gas and star formation inside the inner 2 
kpc of the two simulations. 

What is causing the difference in cold gas in the center? 
There are several possiblities: 

• Cold gas is more efficiently ejected from the early stellar 
feedback simulation 

• Cold gas is replaced by hot gas in the inner 2 kpc. 

• The increased thermal pressure supports the gas at dis- 
tances from the center of the galaxy further than 2 kpc. 
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Figure 12. The evolution of the inner 2 kpc of the stars and 
gas in the simulation for the fiducial (red) and 120% SN Energy 
(cyan) models. The solid lines show the star formation rate inside 
2 kpc throughout the simulation. The star formation rate scale is 
on the left y-axis. The dashed lines show the cold (T< 3 X 10 4 K) 
gas mass evolution within 2 kpc. The dotted lines show the hot 
(T> 3 X 10 4 K) gas mass evolution within 2 kpc. The gas masses 
are given on the right y-axis. 



Rather than the ESF simulation more efficiently ejecting 
gas, a look at the evolution of the galaxies^ shows that the 
120% SN Energy ejects gas at higher velocity and to larger 
radii than the early stellar feedback simulation. The stronger 
gas ejection from the 120% SN Energy simulation is also ap- 
parent in the star formation history. The simulations with- 
out early stellar feedback have less star formation after z ~ 1 
than those with early stellar feedback because they have less 
gas remaining in the central region. 

Figure [12] also shows the mass of hot gas (T > 3 x 10 4 
K, dotted lines) in the central 2 kpc as a function of time. 
The hot gas mass is quite similar for the two simulations, so 
it is not T > 3 x 10 4 gas alone that is supporting gas from 
collapsing into the center. 

Thus, we are left with the possibility that T < 3 x 10 4 
K gas provides enough pressure to support gas to distances 
greater than 2 kpc. To check this possibility, Figure [T3l shows 
the cumulative mass profile of gas at z = 2.25 in the fiducial 
simulation compared with MUGS and the simulation that 
uses 120% SN Energy, z — 2.25 is a time when there is a large 
discrepancy between the mass contained inside 2 kpc for the 
simulations with and without early stellar feedback. There 
are also no satellites within 50 kpc of the disk at z = 2.25. 
Figure [13] shows that there is less gas inside the central 2 
kpc of the fiducial simulation than the 120% SN Energy or 
MUGS simulations. This rapidly changes outside 3 kpc as 
the early stellar feedback simulation has twice as much gas 
there. The increase extends from 3 kpc well beyond r V i r . 
In the case of the 120% SN Energy simulation, the missing 
gas has been ejected out that far as the cumulative mass 
catches up to the fiducial simulation at 300 kpc. In MUGS, 



b http: / /www.mpia-hd.mpg.de/~stinson/magicc/movies/c.ltd.05rp.l/close.mp4 
versus http: / /www. mpia-hd.mpg.de/~stinson/magicc/movies/esnL2/close.mp4 
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Figure 13. The cumulative gas mass profile at z = 2.25 for the 
fiducial simulation compared to the 120% SN energy. The 120% 
SN energy simulation has up to ten times more mass at radii 
inside 3 kpc. In the fiducial simulation, this mass is supported 
outside 3 kpc. 
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Figure 14. A comparison of the distribution of mass in tem- 
perature density phase space within 2 kpc of the galaxy centre. 
The cyan contours represent the 120% SN Energy simulation dis- 
tribution. The red contours show the fiducial simulation. There 
are ten contour levels placed at logarithmic mass intervals. Both 
the temperature axes have been collapsed into one-dimensional 
histograms along each axis. 

however, the cumulative mass well below because the gas 
has been turned into stars. 

To better understand why the fiducial simulation con- 
tains so much less gas inside 2 kpc than the 120% SN En- 
ergy simulations, Figure [14] shows the temperature and den- 
sity distributions for gas within 2 kpc in the two simula- 
tions. Both the density and temperature distributions are 
projected into one-dimensional mass histograms in the top 
and right panels. The temperature mass histogram shows 
that the peak of the gas mass distribution is at ISM condi- 
tions, T < 3 x 10 4 K, in both simulations. It also shows that 



the fiducial simulation has little gas at T < 10 4 K, whereas 
the 120% SN Energy simulation has a significant amount. 
This discrepancy shows how the early stellar feedback is hav- 
ing an effect similar to UV radiation by maintaining gas at 
T > 10 4 K. The UV heated gas at 10 4 K will provide pres- 
sure support whereas the colder gas in the 120% SN Energy 
simulation provides little. Consequently, the 120% SN En- 
ergy simulation contains more high density gas inside 2 kpc 
including gas up to 100 cm -3 . 

The simulations have a similar amount of gas at T > 
3 x 10 4 K, but the hot gas is lower density in the fiducial 
simulation. This shows that the thermal feedback gas is re- 
leased in lower density regions in the fiducial simulation and 
has an easier time adiabatically expanding. 

Our analysis shows the necessity for multiple forms 
of stellar feedback. Adding SN Energy alone simply blows 
stronger winds, but does not create sufficient 10 4 K gas to 
provide pressure support in the disk that keeps low angu- 
lar momentum gas out of the center of the disk. Without 
this feedback, the bulk of star formation happens too early. 
Thus, it seems that winds alone are insufficient to prevent 
the formation of galaxies with high central concentrations. 
It is also necessary for stellar feedback to create a warm 
component of the ISM. That said, the fiducial simulation 
ex hibits the sam e angular momentum distribution described 
in lBrook et alj (2012a), so outflows still play a key feedback 
role in shaping the galaxy. 



5 CONCLUSIONS 

We present a parameter study that varies the strength of 
stellar feedback in simulations of galaxy formation. As a 
constraint, we use the stellar mass-halo mass relationship. 
To reduce star formation enough to match the relationship, 
early stellar feedback from massive stars was required before 
they explode as supernova. 

Our fiducial model best fits the stellar mass-halo mass 
relationship, the evolution of that relationship, has a flat 
rotation curve and an exponential surface brightness pro- 
file with a modest bulge in the center. In contrast to many 
previous cosmological galaxy formation simulations, most of 
the star formation occurs after z = 2. There are two peaks of 
star formation that correspond to the two significant merg- 
ers to the galaxy. Both of these mergers are gas rich. 

We find that the mass of stars formed at z = is very 
sensitive to the amount and timing of stellar feedback em- 
ployed. Models that form too many stars follow a common 
pattern where they turn half of their baryons into stars, 
the other half is hot gas. These galaxi es also form many o f 
their stars at higher redshift than what iMoster et al.l (|2012h 
find when they compare high redshift luminosity functions 
with N-body simulations. Galaxies that form too many stars 
have bright central concentrations of stars. This is reflected 
in galaxy rotation curves as a high central peak. Galaxies 
that form the right amount of stars have exponential surface 
brightness profiles and slowly rising rotation curves. 

A simulation which does not use early stellar feedback, 
but increases the supernova energy to 1.2 x 10 51 erg ends 
up with a similar stellar mass to the fiducial simulation, but 
does not compare as well with other observed galaxy prop- 
erties. The difference is that our implementation of early 
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stellar feedback keeps the gas above 10 4 K and at densities 
below 10 cm -3 , an effect similar to that achieved by UV ion- 
ization. The low density, warm gas keeps the disk extended 
prior to z = 1 and thus keeps gas out of the central 2 kpc 
of the galaxy. Since less gas makes it to the central region, 
fewer stars form before z = 1 and, therefore, this model does 
not result in the massive central concentration that forms 
without early stellar feedback. 

The early stellar feedback produces galaxies that corre- 
spond to observations in a number of ways. Using early stel- 
lar feedback in dwarf galaxies spanning a wide mass range, 
iBrook et al.l (2012b) showed that the simulated galaxies 
match the disk scale length, gas fr actions, and luminosities 
of observed galaxies. IBrook et al.l (2012a) has shown that 
the reason for such agreement is that outflows redistribute 
low angular mo mentum gas from th e centers of galaxies to 
their outskirts. IStinson et alj (|2012l ) showed that the out- 
flows that redistribute the gas also populate the circum- 
galactic medium with an amount of oxyge n that closely 
match es observations of OVI in the CGM. iMaccio et al.l 
l|2012h showed that these outflows can also change the inner 
regions of dark matter density profiles from steeply rising 
cusps into flat cores. In Kannan et al (in prep), we will show 
how the stellar feedback described here performs in a larger 
sample of lower resolution galaxies in a cosmological volume. 
Early stellar feedback helps keeps gas out of the center of 
galaxies, which leads to forming disk galaxies like those that 
we observe. 
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